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Transport in superconducting heterostructures is very successfully described with qnasiclassical 
Green’s functions augmented by microscopically derived boundary conditions. However, so far the 
spin-dependence is in the diffusive approach included only for limiting cases. Here, we derive the fully 
spin-dependent boundary condition completing the Usadel equation and the circuit theory. Both, 
material specific spin-degrees of freedom and spin-dependent interface effects, i.e. spin-mixing and 
polarization of the transmission coefficients are treated exactly. This opens the road to accurately 
describe a completely new class of mesoscopic circuits including materials with strong intrinsic 
magnetic structure. We also discuss several experimentally relevant cases like the tunnel limit, a 
ferromagnetic insulator with arbitrarily strong magnetization and the limit of small spin-mixing. 

PACS numbers: 73.23.-b, 72.10.Bg, 74.25.F-, 72.25.-b 


A central question in quantum transport is to find 
boundary conditions (BC) defining the current across 
contacts of different materials, like e.g. superconducting 
heterostructures possibly including ferromagnetic mate¬ 
rials. A straight forward approach to this is to match the 
wavefunctions at the interface (see e.g. i)- A generaliza¬ 
tion of this approach containing many transport channels 
and all kinds of different barriers can be formulated with 
a scattering matrix Q. However, due to the huge vari¬ 
ety of effects in quantum transport it is in many cases 
more convenient to work with Green’s functions (see e.g. 
i) instead of a wave functions. These Green’s functions 
are described by the Gor’kov equations [^, which are not 
practically solvable in the presence of interfaces and dis¬ 
order. In the qnasiclassical limit (i.e. the Fermi energy 
is much larger than all relevant energy scales) Q it is 
however possible to split the total Green’s into a ballis¬ 
tic part that is continuous at the interface and can thus 
be described within a scattering approach and a discon¬ 
tinuity 0. In this way, BCs for qnasiclassical Green’s 
functions can be formulated. A generalization to this m- 
proach for magnetically active interfaces is given in [7|, 
leading to a pair breaki^ effect of a magnetic insulator 
due to the spin-mixing [8|. Junctions including strongly 
diffusive superconductors are treated in[^. A more de¬ 
tailed overview can be found in 111, where the need 

for a BG the covers superconductor-strong ferromagnet 
interfaces was pointed out. 


As already implied, for many physical systems like 
superconducting heterostructures mainly two limits are 
considered. The clean (ballistic) limit that assumes the 
transport to be described by ballistic trajectories and the 
dirty (diffusive) limit where the mean-free path is the 
smallest length scale in the system and the transport in¬ 
side the leads is dominated by the impurity self-energy 
from Born’s approximation. A BG for stationary prob¬ 
lems that covers spin-dependent effects based on a scat¬ 
tering approach within the clean limit is given in 0. 
The aim of this paper is to provide a fully derived, general 


BG within the framework of qnasiclassical Green’s func¬ 
tions in the case of strongly diffusive materials. The dirty 
limit BG for the spin-independent case is already known 
since long 13| . In this theory (called quantum circuit the¬ 
ory) a matrix-current in Nambu-Keldysh space [mil is 
defined for which generalized Kirchhoff’s rules apply as 
a direct consequence of the discretization of the Usadel 
equation El- Observable currents (like charge, energy 
and spin currents) are then defined from the energy inte¬ 
gral over subtraces of the matrix-current. An extension 
of these BG’s that consider spin-dependent effects to first 
order in the tunneling is provided in 17 -ll. Another ex¬ 
tension that doesn’t consider spin mixing but an interface 
polarization is found in [H. However, a BC that treats 
all the spin-dependent effects exactly is still missing. We 
note that further works consider time-dependent scat¬ 
tering 1211 or spin-dependent effects in specific models 

[ 11130 . 

The derivations in this manuscript are mainly based 
on the ideas presented in [l3| and [1^. Thus the results 
represent a generalizations of the Nazarov boundary con¬ 
dition for spin-dependent problems. 

Throughout the paper we indicate all operators with 
non-trivial structure in more than one subspace by their 
decorations according to table HI We will use cr/r/E 
for the Pauli-matrices in Spin-/Nambu-/Direction (of 
propagation)-space respectively. We work in the Spinor- 
basis 'F'l' = (V’j, ip\, — V’4.5 V’t). Thus, we defined the time- 
reversal operator as R = —ia 2 Ro, where for any A the 
operator i?o yields the complex conjugate {R^A = A*). 

In this basis the transfer matrix M. that transfers 
states ((^■^/«)T = {L/R label the 

left/right side of the interface, e//i labels the elec¬ 
tron/hole component) like has the following 

structure in Nambu space: 


M = 




_0 


( 1 ) 


/M^ 0 

0 M^J V 0 -ia2{MY*i(T2 

The transfer matrix is equivalent to the scattering ma- 
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trix and can be found using the unitarity condition. As 
shown in 18| , the transfer matrix M relating the ballistic 
Green’s functions like = Mg^M'^ is defined as 


M = 


0 

0 


( 2 ) 


Thus assuming the decomposition = (Ad®)' + 

for the electronic transfer matrix, with the 
magnetization unit vector fh, we find the compact form 
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-7-1 
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TABLE I: Explanation of decorations used throughout the 
paper. The symbols / and X indicate that the operator has 
or has no structure in the particular subspace. 


M = M' + M"k. (3) 

find 


Here we defined M' = (7W®)', M" = (Ad®)" and k = 
T3 (g) rha. To find the matrices M' and M", we start with 
the most general form of the electronic transfer matrix 
(see e.g. |l8l|l: 




g»Vn/2 




V’ 


1-T„ -iY„/2 


-i „-iipn /2 


(4) 


The index n labels the transport channel. We as¬ 
sumed the transfer matrix to be already diagonalized 
in channel space 2^ , which is always possible since the 
isotropic Greens functions are structureless in that space. 
Note that the transfer matrix includes channel-dependent 
polarizations of the transmission coefficients P„, spin- 
dependent phase shifts (spin mixing angles) and 

even the spin-quantization axes Explicitly we de¬ 
fined Tn = Tn{l + PnrhnB) and -I- 

((5(^^/^/2)m„(T, while (pn_= and 

The desired form of Ad®, Eq. dS]) is found using the 
matrix structure of fha to rewrite the prefactor of the 
matrix elements in Eq. (|4]) according to 





,_ {t+ — t-ma) 

y/T' 


(5) 

( 6 ) 


with 


T' = T(l-P^) 


(7) 

p+/- = signP°/^^ 

/i(l±Vl-P^) 

( 8 ) 

t+/_ = signP°/i^ 


(9) 


M(i = (M')ii = - ip-s^) (10) 

M" = (Ad")n = -^P^l\-p-c^ + *p+s^) (11) 

= {M'U = (12) 

M" = (Ad")i2 = ^e'>^/^(-t_c^ + tt+s^). (13) 


We defined Cq, = cos(i5q;/4) and Sa = sin(da/4). All other 
components are related by further symmetries (equiva¬ 
lent to the unitarity of the scattering matrix): 

m '2 = {M[,r ; ^22 = «i)* 

M'l = {M[^r ; = (Ad" )*. 

The matrix-current in the isotropization zone on the 
left/right side of the contact is found in O, [13. For 
simplicity we define G = G' = and K = 

K^lR^ where the indices L/R label the left/right side of 
the interface respectively and G^/^ are the quasiclassical 
Green’s functions of either sides. With these definitions 
and defining the conductance quantum Gq = jh^ the 
matrix-current (/ = I^/^) is 

/ = 4GQTr„,,[(l + GK)- 1 (GE 3 ± 1) T ^1, (14) 

with 

id = (Adt)±i(/./^±i^ (^5) 

The index s labels the trace in direction subspace. 
For later usage we mention that due to = 

= E 3 and the normalization condition (G)^ = 1 , 

it is: 

(E 3 X)" = 1 . (16) 

The inverse of the 2x2 matrix in the direction subspace 
in Eq. m is found in terms of components Kij of K as 
follows. We define 


Note that all these quantities (transmission, magnetiza¬ 
tion, polarization, phases and phase shifts) are diagonal 
but non-trivial matrices in the channel subspace. We ex¬ 
pand Eq. ([3]) with the help of Eqs. and m and 


L»-i = (l-bGid)-i 

1 + GPii Gki2 {D-^)i2\ (17) 

Gid21 1 + G^22y V(P-')21 (^-')22; 









































Performing the trace over the s-subspace in Eq. da we 
find: 

/ = 4GQTr„[(D-i)n(G ± 1 ) - ( 6 - 1 ) 22(6 T 1 ) T !]■ (18) 
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using the relation [A 1 , G] 
show that 281 


0 it is straight forward to 


The entries of the 2x2 matrix D ^ are found by direct 
inversion: 

(D-i)ii = (1 + GKii - Gki 2 {l + G^22)“^Gif2i)“^ ( 19 ) 
(D-1)22 = (1 + Gk22 - Gk2i{l + G^ii)-iGi^i2)-i.(20) 

Using the matrix structure of k (Eq. (fTHl) '). i.e. 
kkk 22 = kiik;^^ and kk + ki 2 = kiik;^^^k 22 , 
Eqs. (IT^ and (12(1 can be decomposed according to: 

(6-i)n=Gi-i (21) 

(6-1)22 = A-\k^^^G + G^n^^-iiG) (22) 

with 

G = (6621)-! (1 + G622) = k^G + 6116211 (23) 

A = 62 IIG + G 61162 IIG + 62 I 1622 
+ G 61162 II 622 - G6i2)-1 
= G + GGG. (24) 

Plugging Eqs. (HJ), ([11]), ([13]), and (HH) into Eq. (IT31) . 
replacing the last summand =pl in Eq. (IT8)) by A-^A and 


/ = 4GQTr„[(l±6ii)(62i)-i, 

{(G + 6ii)(62i)-i,G}-1(G±1)]. (25) 

This is the most general form of the desired boundary 
condition, and thus the main result of the paper. To solve 
physical problems, we still need to find the components 
of 6 to be plugged into Eq. (EH). For this we plug Eq. ® 
into Eq. (fT51) (note that detM = 1 ): 

K=((M')^k\M'k^G' + ((M")^k^(M"k^KG'K 
+ ((M')^k\M"k^G'K+((M")^k^(M'k^kG'. (26) 

The last two summands are related by the symmetry 
(((M')t)±i(M")±i)t = ((M")l)±i(M')±i. Plugging 
Eqs. (fTUl) . (fTTl) . (fril) and (fT3]) into Eq. [35] using all the 
symmetries one finds the following for the desired com¬ 
ponents of 6 : 

611 = i^lG' -b K'l^kG'k + K'l[G'k + K'{'*kG' (27) 
621 = X 21 G" -b K”^kG'k + K”[{G', k} (28) 

with 




1 1 


i± 


- {\/l-P2 cos(( 5 (^/ 2 ) -b 6(1-T02_p2 cos(dx/ 2 )} 
+ ^ { 61 - P2 sin(d(^/2) - v^(l - T')2 _ p2 sin (^x/2)} 


(K'^ir = ^ [-P 

{^21 k = [(p+t+ +P-t-)cosiS<j)^/2) ± (p+t+-p_t-)cos(54>^/2) 


-i {{p-t+ +p+t_)sm(S(l)^/2) =F (p+t- - p-t+)sm{S(l)^/2)}] 


{K'kk = [(p+t_ -bp_4)cos((5<))^/2) - i{p+t+ + p_t_)s\n{54)^/2)] 


(29) 

(30) 

(31) 

(32) 


Here we partly used Eqs. m and ®. Due to the structure of the transfer-matrix the right-side coefficients are given by 
^ )^ = Po{Kab ^ Now we Can express the full boundary condition with (0 = k = = 

5(j)' = (j)^/^) 


611 = (G" -b kG'k) 


l-Pk 


T' 


-b — (G' - kG'k) 61 -P2 g±*«(5<b+50')/2 ^(1 _ r')2 _ p2 g±*ft(50-50')/2 


K2^ = - 


I (G' -b kG'k) [p+t+ +p-t-- k{p+t- + p-t+)] 

+ (G' — kG'k) [(p+t+ — p-t-) cos((56/2) ± i(p+t- — p-t+) sin(d(()'/2)]| . 


(33) 


(34) 


This completes the central result of our work. Plugging Eqs. (l33)) and (l34ll into Eq. (1251) yields the dirty limit 
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matrix current through a contact where both sides (G^ 
and G^) can have a non-trivial structure in Nambu and 
spin space and the spin-dependent interface effects i. e. 
spin mixing and spin filtering are treated exactly. It com¬ 
pletes the discrete version of the Usadel equation [1^ to 
yield a complete scheme to describe a large variety of 
superconductor-ferromagnet heterostructures, with ex¬ 
ternal fields under equilibrium and non-equilibrium con¬ 
ditions. 

In the remainder of the paper we discuss several prac¬ 
tically relevant simplifications of the general boundary 
conditions. At first we note that the Nazarov boundary 
condition 0 is reproduced for 6(j) = 6(j)' = P = 0, where 
Kn = (2 - T)G'/T and K 21 = -2y/t^G'/T. Thus, 


with [G", ( 2 ( 1 ^ 2 ) reduces to 

the well-known result 


^ = 4Gq^ 

n 


Tn[G',G] 

4 + r„({G',G}-2)' 


(35) 


In many cases tunnel barriers (with or without a mag¬ 
netic structure) are an experimentally realized and we 
will at first take the tunnel limit (T„ <C I) by linear ex¬ 
pansion of Eqs. (l3^ and (|34l) : 

T'Kn/2i«^^iV2i+^n/2i?^ (36) 


with 


= (1 - Pk){G' + kG'k) + (G" - kG'k)VT^^cos {54)'/2) 
k\^ = _i [(1 _ p-^)(G' p kG'k) + \/l-P2(G" - kG'k) 5-/-')/2 


2 1 

^0 _ _ ±i<p 
i\ 2 i — e 


(1 - Pk){G' + kG"k)e =^“^'^/2 ^i-p2(^G' - kG'k) cos (54)'/2) 






(1 - P^)iG' + kG'k)e^''^^‘^l'^ P \J\-P‘^{& - kG'k){cos {54)'12) p iPsuv{5(4'12.)) 


(37) 

(38) 

(39) 

(40) 


Plugging Eq. (|55)) into Eq. (1^ results in 

/ = 4GQTr4K?i(.A°i)-\ {^n(^2i)-\ G}-'(1 - {(G(l - P^) P kl, - k^^,)-^kl,){k^^^)-\G}x 
x{kU^i)-\G}-^T){i ±G)] + [(±(1 - p2) + kl, -kU^^)-^kl,){^^)-\ 
{kUK^^i)-\G}-\l±G)]T. (41) 


Using the relation Va G C : e *''“[G",k] = [G",k]e*'‘“ 
one finds 

^ii(-^2i)”^ = (42) 

kl^ P = 

= ^^Jl-p■2[G', k](I - Pk) sin (5072) = ±k^, (43) 

and thus 

/ = 4Gq G}-7G ± I)] 

n 

+T„([e^“" G}-7(G(1 - P^) ± k^) 

x(IG27)-i, G}{e^*"" G}-1(G ± 1 )] 

-[(1 - Pi P kl){kl^)-\{e^^^- G}-i(G ± 1)]).(44) 

Eq. (l44ll together with Eqs. (1431) and (l3^ represents the 
tunnel limit of the full boundary condition (eq. (|25l) l. 
In the case of a ferromagnetic insulator = 0 only the 
first line in the Eq. (1441) remains. Further simplifications 


are found assuming small spin-mixing. We note that 
({G', k}{k -P)p [G', k]kVl - P2) = 

{G", k}{k -b P) + [G', k]k^l — P^ 

4(1 _P2) • Z t>l 

By expansion of Eq. (I44p to first order in the spin- 
dependent phase shifts (50„, 50)^), using Eq. (HSl) we 
find: 

2/ = Gq ^ [±r„ ((I + ^r^)G' + P„{G", kn} 

n 

+ (1 - V'l - PDknG'kn) - 2l 54)n, G . 

(46) 

This is the boundary condition as it is used in [l^, 
[ 29 II to calculate thermoelectric effects in ferromagnet- 
superconductor heterostructures. 

In conclusion we have found a new boundary condi¬ 
tion for isotropic quasicalcassical Greens functions that 
generalizes the Nazarov boundary condition in the quan¬ 
tum circuit theory to spin-dependent problems. It is 
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presented in Eq. (ESI) supplemented by Eqs. El) and 
(1341) . The new boundary condition captures e.g. non- 
collinear ferromagnetes and triplet superconductors as 
terminals and fully accounts for arbitrary transmissions, 
polarizations and spin-mixing effects of the scattering re¬ 
gion. In the tunnel limit the boundary condition can 
be simplified to Eq. (|4^ . still fully accounting for arbi¬ 
trary spin-mixing. This case also comprises ferromag¬ 
netic insulators by setting T„ = 0 given by the first line 
of Eq. (l4^ . The next drastic yet non-trivial simplifica¬ 
tion of the boundary condition is to linearize in the spin- 
dependent phase shifts given in Eq. (|46)) . In this case 
the boundary condition takes a form of a commutator 
(similar to [l^ and in [l^ ) and the transport coefficients 
just enter in averaged interface parameters T„, P„ and 
S4>n^^- Our results pave the way to investigate spin- 
dependent charge and heat transport and spin transport 
in superconducting heterostructures as well as unconven¬ 
tional correlations like odd-frequency triplet pairing us¬ 
ing the most general boundary condition based on a fully 
microscopic derivation. Another application is to verify 
the details of microscopic models describing the intrin¬ 
sic magnetic structure and magnetization behavior of e.g 
ferromagnetic insulators. 
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